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Abstract 

n 

The reason behind the discrepancy between the phase diagrams of our earlier work |l( and 
the comment of Sheehy and Radzihovsky 2] is discussed. We show that, in contrast to what is 
claimed in j^J, the requirement of positive susceptibility is sufficient to rule out states that are local 
maximum of the free energy (as a function of the order parameter A). 

PACS numbers: 03.75.Ss, 05.30.Fk, 34.90. +q 
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It is now widely accepted that, for an attractive s-wave interaction, an equally populated 
Fermi gas at zero temperature smoothly crosses over from the BCS to the BEC regime. In 
our recent paper jlj], we considered a two component Fermi gas with unequal populations 
under a wide Feshbach resonance. We showed that the uniform state must become unstable 
at some intermediate coupling strength since it has either negative superffuid density or a 
negative susceptibility. Therefore we demonstrated that the smooth crossover known for the 
equal population case is destroyed (independent of the ansatz what actually replaces the 
unstable states). A phase diagram was then constructed by indicating where the uniform 
state was found to be unstable (reproduced in Fig 1 here as the region between the dotted 
lines) . 

In a recent preprint, Sheehy and Radzihovsky |2| reinvestigated this phase diagram by 
a different method, and found discrepancy with our earlier results (see also jj, 0, |f|) 
In particular, they found that the instability region occupies a larger area than ours on the 
BEC side of the the phase diagram. They suggest that the susceptibility criterion [j] we 
used is not sufficient to rule out unstable states that are actually relative maximum of the 
free energy Q with respect to the pairing potential A. 

We here would like to clarify the reason causing the discrepancy between the above two 
works . We now believe that our phase diagram in [l| is incorrect. However, the reason 
was in fact due to the inaccuracies of the numerical method we used there. Further, we 
argue that the susceptibility criterion is able to rule out states that corresponds to relative 
maxima of fl (such as those depicted in Fig 2 of_Q), in contrast to what is claimed in Q]. 

First we comment on the method we used in |l[ . A stable state must have all eigenvalues 
of the susceptibility matrix dn a /dfi a i positive (here a and a' =j, \ for the two species, and 



A*T,(I) = fJ- ± h are their chemical potentials and n a are the densities). In 1], we solved, for 
the uniform states, the dimensionless average chemical potential fi = fi/tF and potential 
difference h = h/ep as a function of the dimensionless population difference fid = n^/n and 
coupling constant g = 1/kpa. (n = n^ + n^ rid = —n\). The inverse susceptibility matrix 
d[i a jdn a i can be expressed in terms of the above obtained functions fZj. It can further 
be shown [t| that a necessary condition that follows is that [dhjdh^j must be positive 
for stability, i.e., that the slope of h versus fid must be positive at fixed scaled coupling 
constant g. Since the functions fl(n d .g), h(fi d . g) were already available (sec. e.g.. Fig 1 
of [J), we chose to evaluate this matrix and its eigenvalues numerically from these data, 
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and constructed the phase diagram in Unfortunately, as we only found out later, the 
numerical accuracies required to carry out this method is very high, and we erroneously 
concluded in |l| that one of the eigenvalue changed sign at the same position as where 
dh/dfid changed sign (footnote [13] of We have now re-evaluated the positions where 
(one of) the eigenvalue of this matrix changes sign more accurately, and the result is as 
shown a 8 large circles in Fig 1 in this reply. As far as we can tell, this cnterion yields the 
same line (dot-dashed) as where d 2 Vt/dA 2 changes sign 8]. Thus on the left of this line, 
we find that our stable state has a positive susceptibility matrix as well as corresponds to 
a relative minimum of the free energy Q, whereas on the right the unstable state identified 
has a susceptibility matrix with negative eigenvalue as well as corresponds to a relative 
maximum of the free energy Q. Therefore, the susceptibility criterion is able to identify 
states that are relative maximum o 
in contrast to what was claimed in 

The most transparent way to see the last statement is to consider the expression 9. If 
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the free energy and concludes that they are unstable, 
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Note that the matrix (^§jj 3 ^j has eigenvalues that are positive and finite (e.g. pf), and 
the matrix (^Sr) |-5aH has non-negative eigenvalues. Consider now a stable state 
on the BEC side at large g. This state has positive susceptibilities and is a free energy 
minimum, with d 2 Q/dA 2 > 0. Consider now decreasing g, thus moving towards resonance. 
At the point where this state is no longer a free energy relative minimum but just becomes 
a relative maximum, d 2 Q/dA 2 < and small. Thus the second term in eq and there- 
fore the susceptibility matrix itself, necessarily has a large and negative eigenvalue at this 
point. Therefore, as we claimed above, the susceptibility criterion is always able to identify 
the situation where, as coupling constant changes, a relative minimum becomes a relative 
maximum, (see also 0,0]) 

We agree however with Sheehy and Radzihosky Q] that our susceptibility criterion may 
not always protect us from instability caused by a first order transition. But this can oc- 
cur only when a new free energy minima arises at some other A, rather than having the 
original minimum turning into a relative maximum. In general, our criterion can fail when 
one can come up with a better ansatz for the state (at the same n and n d ) than the one 



originally investigated, or there is a physical instability not contained in the susceptib ility 
being evaluated. This happens in particular on the BCS side of the phase diagram [11]. 
To the right of the dotted line on the BCS side, the normal state is a free energy relative 
minimum. However, as suggested by the negative susceptibility matrix and superffuid den- 
sity (for the latter, except near « 1) of the uniform state to the left of this dotted line 
[l], a phase separated state and a state with finite pairing momentum (FFLO state) are 
potentially more stable than the normal state in a region near and to the right of this dotted 
line. One can straightforwardly evaluate these instability lines corresponding to these two 
ansatz. For phase separation, one finds the point where the free energy of the normal state 
becomes higher (when g is increased, i.e., moving from right to left in Fig 1) than that of the 
completely paired superfluid state (at the same chemical potentials fi and h of the normal 
state) We evaluated the latter free energy by an integration over coupling constant Q]. The 
FFLO instability line of the normal state can be found by solving the Cooper problem at 
finite pair momentum q, again at the chemical potentials of the corresponding normal state. 
These lines are also shown in Fig 1 here Q|. The diag ram constructed from these transition 
lines then agree with [J| (see however the remark |l l| ) . 

In conclusion, we have resolved the discrepancy between the phase diagrams in |l| and 
. The reason for the earlier disagreement was clarified. 
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While we were finishing the present manuscript, we noticed another preprint 14| dis- 
cussing the same topic as here. 
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2112-M-001-002 (SKY). 
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we have used the fact that in three dimensions, €f oc n 2 / 3 and kp oc ra 1 / 3 . The condition that 



all eigenvalues of B aa i are positive requires that its expectation value with respect to the 
vector (1, —1) be positive, and hence B^ + B^ — B^ — B^ > 0. Using the above expression 



for B aa i , we get the necessary condition ( > 0. 
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where E^s are the quasiparticle 
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FIG. 1: (color online) The phase diagram, with the region where the uniform state being unstable 
shaded. The stable region to the BCS (right) side corresponds to the normal state, whereas the one 
on the BEC (left) side corresponds to a gapless superfluid. The boundaries for the unstable region 
reported in [1] are indicated by the dotted lines. On the BEC side, the dot-dashed line represents 
where d 2 Q/dA 2 changes sign. Circles are where an eigenvalue of the susceptibility matrix changes 
sign. On the BCS side, the dashed (blue) line represents the location where the normal state 
becomes unstable towards phase separation. The full (red) line represents the finite q Cooperon 
instability of the normal state. 
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